Legitimacy of wave-function expansion 
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Abstract: In this letter we investigate the common procedure in which any 
wave function is expanded into a series of eigenfunctions. It is shown that 
as far as dynamical systems arc concerned the expanding procedure involves 
various mathematical and physical difficulties. With or without introducing 
phase factors, such expansions do not represent dynamical wave functions. 
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A common concept, which is originally due to Dirac and Neumann and 
has been accepted as one of the essential principles of quantum mechanics, 
is that any wave function can be represented by a vector in the linear infinite- 
dimensional vector space, the Hilbert space, whose basis is formed by a complete 
set of orthogonal eigenfunctions [1-4]; namely a wave function can be expanded 
into 

£0,(i)¥„(r), (1) 

n 

where n stands for a set of quantum numbers labeling the eigenfunction and 
C n (t) a pure time-dependent complex number. The expanding procedure is 
also called the principle of superposition in some textbooks [4]. Adopting the 
expanding procedure, we can write a dynamical, nonstationary, wave function 
in the form (with no spin regarded) 

£C n (t)exp(-i^t)* n (r), (2) 

n 

where e„ is the eigenenergy of the system before a certain initial time to and 
ty n (r) the eigenfunction in terms of the Hamiltonian H(to)- By substituting (||) 
into the time-dependent Schrodinger equation and making use of the orthogo- 
nality of the eigenfunctions, we arrive at a set of coupled ordinary-differential 
equations 

ih—rr = X! Ci(t)(Hi) n i exp(iuj n it), (3) 
i 

where Hi = H(t) — H(t ) represents the variation of the Hamiltonian (not 
necessarily small) . 

It is almost unanimously believed that (|l|) or (|J) formally represents the 
real wave function, and that the set of ordinary-differential equations (||) is 
equivalent to the corresponding Schrodinger equation. Furthermore, the value 
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of \C n \ 2 is, in Dirac's perturbation theory, interpreted as a transition probability 
from one quantum state to another. 

While the expanding procedure outlined above seems stringent and flaw- 
less, many difficulties that can be traced back to it have constantly bothered 
people in the community. Several decades ago it is noticed that the resultant 
formula derived from the expanding procedure in Dirac's perturbation theory is 
not gauge-invariant [5] . Later, some in the community found that coefficients of 
such expansion do not represent transition probabilities [6]. To cope with these 
difficulties, people kept on coming up with remedies. The preferential gauge, in 
which the vector potential vanishes whenever the electromagnetic field becomes 
zero, was proposed after the debate on the gauge aspect of Dirac's perturbation 
theory [7]. Furthermore, it was suggested that a certain phase factor should be 
introduced to the wave function before using the standard expanding proce- 
dure [8]. 

In our view, those remedies have obscured, to some extent, the essence of 
the problem[9]. Actually, we can easily see the situation is unsatisfactory with 
or without these remedies: if we adopt the expanding procedure as it is, we take 
the risk to confront with other physical principles such as the gauge invariance; 
if we believe that introducing extra additives can make the procedure work as 
it should, the expanding procedure becomes more a recipe than a basic physical 
principle. 

The objective of this letter is to directly investigate the validity of the ex- 
panding procedure. It will be illustrated that as far as dynamical systems are 
concerned the expansion form expressed by (|I]) or (^) does not represent a true 
dynamical wave function. If we forcefully use such expansions, abnormal things 
will emerge. Furthermore, it will be shown that introducing extra phase factors 
does not really help. 

It is very obvious that conclusions concerning this subject are crucially im- 
portant since many theories and many discussions in the literature are explicitly 
or implicitly related to the expansion principle discussed here. 

Before starting with our main discussion, it is appropriate to recall how the 
expanding procedure works for a stationary system. For such system, we have 
a time-independent Hamiltonian Hq and eigenfunctions \l/ n (r), and any wave 
function of the system can indeed be expressed by 

*(i,r)=^C„exp(-^t)* n (r), (4) 

n 

where C^'s are pure numbers independent of time and space. Note that the 
normalization condition 

£ic»i 2 = i ( 5 ) 

n 

holds for (||) and it ensures the unity of the total probability and the convergence 
of the expansion. 
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We now try to determine whether or not a wave function of nonstation- 
ary system can be similarly expanded. We will first investigate mathematical 
arguments concerning this issue. 

Consider a quantum system that has a stationary Hamiltonian Hi before the 
initial time U and another stationary Hamiltonian Hf after the final time tt. 
During the time period ti < t < tt, the system involves a dynamical change. In 
using an expansion of the form (yj) to describe the wave function of the system, 
we need to select a set of eigenfunctions we wish to use. We may choose the set 
of eigenfunctions of Hi which we denote by or we may equally choose 

the set of eigenfunctions of Hf which we denote by { x f' n/ }. Referring to the 
discussion on the stationary system, we know that the wave function can be 
expressed as, before U, 



exp 



("ttO 



(G) 



and after tf 



(7) 



In harmony with the expansion (||), we may simply use the eigenfunctions {^n;} 
to expand the wave function throughout the dynamical process. To ensure all 
the things go smoothly, we need to verify that any one of {^nr } can be expressed 
by {*& ni }. If this is possible, we will say the two sets of eigenfunctions, or the 
two states, are compatible. 

It is found that for some dynamical systems the final state is incompatible 
with its initial state. As an explicit example, we take a look at a particle that 
is trapped by a magnetic field initially and then moves in a free space since 
the magnetic field disappears. An eigenfunction for the particle in a uniform 
magnetic field is [4] 



/27T 



R n ,i{p)e 



il<p ik x z 



(8) 



where 



with 



RuAp) = T{ ^llf l e(-&F(-n, \l\ + l,p 2 /2« 2 ) 



(\l\+n)\ 
2\ l \n\ 



and F is the confluent hypergeometric function 

„, , a u a(a + 1) u 2 

F(a, 7 ,u) = 1 + -- + - 
7 1! 7(7 + 1) 2! 



(9) 
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On the other hand, an eigenfunction for a free particle motion is [3] 

(87r 3 )- 1/2 exp(ik-r). (10) 
If we particularly choose a plane wave as 

* fe ,(r) = C"exp(ifc^), (11) 

and expand it into 

* fc /(r) = J2 ^,1,*. (12) 

n,l : k z 

where ^ n ,i,k z is an eigenfunction of the type (||). The coefficient C nt i^ z has three 
parameters. For the purpose of this letter, we set I — and will be concerned 
with the convergence of the series in terms of the index n only. Along this line, 
we obtain 

C n , .k z oc S ( k *- k 'J J e~&F{-n, 1, P 2 /2al)pdp. (13) 

It happens that the expression above can be evaluated analytically and the 
result is, with the delta-function disregarded, 

Cn+i,o,fc a = C ni o,k z ■ (14) 

which means that the series of the form (jl^) is not convergent in terms of the 
non-negative integer n. 

Actually, the difficulty can be seen just after we write down the two sets of 



eigenfunctions, (|10|) and (J8|) . In the x — y plane, eigenfunctions ( J10| ) take finite 
values at the infinity while eigenfunctions (||) tend to zero at the infinity. It is 
obvious that these two kinds of eigenfunctions cannot express each other. Also 
note that the example presented here possesses a certain degree of generality. If 
we deal with a situation in which a magnetic-type field is applied to an atomlikc 
system (having discrete and continuous spectrum), we will encounter similar 
difficulties. 

If one thinks that the argument presented above is not general enough, look 
at the following argument in which we show that the normalization of wave 
function is broken down with the expanding procedure employed. 

We first write down the formal integration of differential equations (|J). Let 
the disturbance time interval be — to, where to and ijv denote the initial time 
and the final time at which the perturbation turns on and off respectively. Slice 
the time interval into N small equal increments 

to, tx, t2, , £jv, (15) 

where — t{ = (t^ — to)/N = At. The coefficients Ck can be evaluated by 

Cfc(ti) w C k {t ) - (i/?i)ECfc'(*o)(-ffi)fcfe'e i ^' to Af 
k' 

C k (t i+1 ) w C k (U) - {i/h)'£C k ,(ti)(H 1 ) kk ,e iu »-' t *At (16) 

k' 
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According to Euler and Cauchy[10], if we slice the entire time into a sufficiently 
large number of intervals or if we are just concerned with the system's behavior 
within a sufficiently short time, the integral expression ( |l6| ) can be accepted as 
an accurate one. 

Suppose that the system is initially in one eigenstate, namely we assume 

C s (t ) = l, C k (t o ) = 0(k^s), (17) 
then we find that at the later time t\ 

C a (t 1 ) = l- l -{H 1 ) ss M. (18) 

Since H\ = H(t) — H(to) is a self-adjoint operator, (Hi) ss is a real number. 
(For a uniform electric field, we can assume H\ oc x and x ss is just the average 
x-position in the initial state.) Then, (|l8| ) tells us that 

\C s (h)\ 2 > 1. (19) 

For other coefficients, we have 

J2\C k (h)\ 2 >0. (20) 

k^s 

Eqs. @ and @ lead us to 

^|C fe (tx)| 2 >l. (21) 

k 

The equality sign holds only when all the matrix elements of Hi are zero and 
the system gets no disturbed. In that case, we investigate the value of the same 
expression at the time ti- In such a way, we will surely find 

J2\C k (U)\ 2 > 1. (22) 

k 

This clearly shows that the expansion of the form (||) or (|l|) is not the correct 
formal solution of the time-dependent Schrodinger equation. (In this letter, we 
take it for granted that the Schrodinger equation preserves the unity of the total 
probability.) 

One interesting question arises immediately: what does a dynamical wave 
function look like? The discussion above has suggested that such wave function 
must involve coupled timc-and-spacc dependence, which cannot be handled with 
the variable-separation technique. (Finally, incompatible functions may also get 
involved.) Functions that can be expressed by ([[]) or (|J) constitutes a special 
type of functions; but, paradoxically, none of them represents a dynamical wave 
function. 
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We now turn to physical arguments of the issue and in passing we will study 
effects of introducing phase factors. 

Before doing that it is enlightening to recall the well-known concept in quan- 
tum mechanics that the gauge fields, Hamiltonian of system, phase factor of 
wave function and operators representing observables are not uniquely defined 
and they covary in accordance to (m = Q = c = 1 in this letter) 

(A, $) — > (A' = A + V/, $' = $ - d t f) 

*(t,r) -> e W' r >¥(t,r) y ' 

L(p-A,r)->L(p-A',r), 

where f(t,r) is an arbitrary differentiable function of time and space. If and 
only if the covariance relationship given above is observed the physical outcomes 
of the formalism will be strictly gauge-invariant. 

We can immediately see that our expansion principle will suffer in this re- 
spect. The third equation of (J23) illustrates that even if the wave function of 
interest could be expressed by (Q) or ^ in a special gauge, the wave function 
in a general gauge would involve an extra nontrivial phase factor e l ^ l ' r \ Such 
nontrivial phase factors pose a serious difficulty for the standard expanding pro- 
cedure expressed by (|l|), (|j) as well as the differential-equation set (||). This 
partly explains why the time-dependent perturbation theory encounters gauge 
difficulties. 

At this point, it is in order to point out that the standard expansion proce- 
dure is in conflict with the covariance relationship expressed by (|23|). Consider 
a quantum particle that is bounded in a mechanical well. Its state before t = 
is assumed to be 

#(t<0)=exp(£^tj#„(r). (24) 

At and after t — 0, an electromagnetic disturbance is applied to the system. If 
we admit the standard expansion principle, the wave function after t = must 
be in the form 

*(t>0) = ^a(*)exp(i|t)*j(r), (25) 
i 

where 

C a (0) = l, Q(Q) =0 (fcrijig). (26) 

We note that the quantum system has its own inertia (the heavier the particle 
the larger the inertia) so that all coefficients in the expansion, as well as physical 
observables associated with the system, must vary tardily in comparison with 
the variation of the disturbance field. The velocity operator reads according to 
the standard definition 

v = p A = -ffiV - A. (27) 
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where A is the vector potential chosen to represent the electromagnetic distur- 
bance. Before t = 0, the average velocity of the particle is 

v = (¥ a |-tftV|*.). (28) 

At the instant t = + e, where e is a truly small time interval, the average 
velocity becomes 

v=(* 5 |-iW|* s )-(* s |A|* s ). (29) 

The trouble with the expansion (^5j) is now obvious: not only that the aver- 
age value of the velocity is gauge-dependent but also that its value varies as 
promptly as the vector potential. Actually, all observables, including the tran- 
sition probability, suffer from the same difficulty with the use of the expansion 

There remains one final question concerning whether or not the expression 

e i/M £a,(t)$ )l (r) (30) 

represents a true dynamical wave function. Our answer to it is a negative one. 
If this expression represented a dynamical wave function, we could find a special 
gauge under which 

£<7 n (t)* n (r) (31) 

would represent a dynamical wave function. According to the arguments pre- 
sented in this letter, particularly those independent of gauge, this is not at all 
possible. 

In conclusion, it has been shown that the common expanding procedure, 
or the so-called principle of superposition, holds only for stationary systems. 
However, we notice that in the real life many quantum theories concerning 
nonstationary dynamical processes have assumed the general validity of it [11]. 

Discussion with Professors R. G. Littlejohn and Dongsheng Guo is gratefully 
acknowledged. This work is partly supported by the fund provided by Education 
Ministry, PRC. 
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